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Abstract. We argue that the natural way to generahse a tensor network variational class to 
a continuous quantum system is to use the Feynman path integral to implement a continuous 
tensor contraction. This approach is illustrated for the case of a recently introduced class of 
quantum field states known as continuous matrix-product states (cMPS). As an example of the 
utility of the path-integral representation we argue that the state of a dynamically evolving 
quantum field admits a natural representation as a cMPS. An argument that all states in 
Fock space admit a cMPS representation when the number of variational parameters tends 
to infinity is also provided. 



1. Introduction 

It is a remarkable fact that a generic state of a quantum system of moderate dimension is 
unlikely to be physically realisable. This is because physically accessible states occupy only a 
tiny submanifold of the full hilbert space H, as follows from a simple counting argument (see 
e.g. [I3]). The submanifold of physically accessible states has been christened the physical 
corner of Hilbert space; states in this corner exhibit highly nongeneric features: a quantum 
state chosen at random from 7i is massively entangled and locally featureless [4J while, in stark 
contrast, the manifold of physical states realized in Nature exhibit short-range correlations, 
and often obey entropy area laws [5]. 

The fact that physically realisable states are not generic has important practical conse- 
quences for the calculation of observable properties. Instead of parametrising a physical state 
using a complete (exponentially large basis) of hilbert space, entropy area laws suggest that 
there is a more efficient description where only the physical degrees of freedom participat- 
ing in short-ranged correlations are used. Such a parametrisation, if found, would allow 
unprecedented insights when combined with the variational method. 

A compelling example of an efficient parametrisation of the physical corner is found in 
the study of one-dimensional strongly interacting quantum lattice systems [HI [71 [T]. The 
archetypal variational class here is the set of matrix product states (MPS) [6]. The variational 
method applied to MPS, synonymous with the density matrix renormalisation group [231 IT6l 
fT5] . has enjoyed unparalleled success in the study of the equilibrium and nonequilibrium 
dynamics of quantum lattice systems. 

The MPS class has now been extended far beyond the one-dimensional setting to arbitrary 
quantum lattice systems where they are now understood as examples of tensor networks. The 
two most notable tensor networks are the projected entangled-pair states (PEPS) [17j and the 
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multiscale entanglement renormalisation ansatz (MERA) 120] . Both of these classes have 
been apphed with considerable success to the study of strongly correlated physics in two 
dimensions and beyond for noncritical and critical systems. Tensor networks such as PEPS 
and MERA are not yet able to explain all the physics of low-dimensional lattice systems, but 
already many new insights have been obtained. 

While the MPS and MERA classes are naturally defined in the presence of a lattice, they 
have both [TSl [T21 El E] been generalised in a natural way to continuous quantum systems. 
Here a basic challenge is that a discrete collection of degrees of freedom is replaced by a 
continuous collection. This means that, for a tensor network state, a contraction of a discrete 
collection of tensors needs to be replaced by the "contraction" of a now continuous set of 
tensors. Here we argue that the most natural systematic way to achieve this is to replace 
the tensor contraction with a path integral over some now continuous auxiliary degrees of 
freedom. Thus the core objective of this paper is to obtain and study such a path-integral 
representation for the cMPS class. This representation will then be used as the basis for 
higher- dimensional generalisations in a forthcoming paper. 

2. Background: matrix product states, tensor networks, and path integrals 

Here we review the MPS class and sketch some of its properties. Our intent is to make this 
paper accessible to those with a diversity of backgrounds, so we have tried to provide all the 
necessary prerequisite material and references needed to follow our argument here. Readers 
with a familiarity with MPS and the DMRG are invited to skim this section lightly to fix 
notation. 

We begin by recalling that any bipartite pure quantum state iV'^g) admits a Schmidt decom- 
position IV'^b) = Ylik V%;|^)^® 1^)23 for some set of local bases of A and B. For pure states IV') 
of one-dimensional quantum spin systems A1A2 ■ ■ ■ An with local dimension d, we may per- 
form a Schmidt decomposition iteratively on the bipartitions [^1,^'^], [(^1^2)5 (-^1-42)'] , ■ ■ ■ 
[A'n, An] (where X' is the complement of X and [X, Y] denotes the particular bipartite split) 
to obtain the MPS representation [21] 

d-l 

(2.1) 1^^)= Yl MA^'A^'---A^''M\jiJ2---Jn). 

Here A-"", jyt = 0, 1, . . . , c? — 1, is a collection of d matrices of size -D^-i x D^, {ujl\ is a row 
vector of dimension Dq, and is a column vector of dimension Dn- The dimensions Dj. are 
called the bond dimensions of the MPS and characterize the degree of entanglement entropy 
across a cut at site k. This construction shows that MPS are an expressive class, meaning 
that any state may be represented as an MPS with a sufficiently large choice of the D^s (the 
argument applies to any pure state). However, in most implementations we simply assume 
that the bond dimension is constant and truncate it at some value = D, which acts as a 
refinement parameter for this class. 



Matrix product state representations (2.1 ) possess several remarkable properties. The first, 
and most important, is that they provide an efficient parametrisation of naturally occurring 
states [ini [m [13 US]; it is now understood that MPS can very efficiently represent both 
the ground states of models with a spectral gap and also the non-equilibrium dynamics of 
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any quantum spin chain. The second property is that they possess an entropy area law [5], 
meaning that the von Neumann entropy of any contiguous block of spins is bounded above 
by a constant, i.e., the size of the boundary. Another important property of MPS is a gauge 
degree of freedom, so they supply an over-complete parametrisation. 



A matrix product state (2.1) is an example of a quantum state known as a tensor network 
state (TNS). To define a TNS one first associates a finite graph G = {V, E) with the quantum 
system where the physical degrees of freedom, which are quantum spins of dimension ci, live 
on the vertices V , and the edges E encode auxiliary degrees of freedom. To each vertex v 
we associate a tensor A^^ „ „ with + 1 indices, where is the degree of the vertex v. 
Each index ae is associated with a corresponding edge e & E incident with the vertex v and 
takes values 1,2, .. . ,De; these are the auxiliary bond indices. The index is the physical 
index and takes values 0,1, . . . ,d — l. The TNS corresponding to this arrangement of tensors 
is then given by 

(2.2) 1^) = 



where C denotes the contraction of all the auxiliary indices, i.e., each pair of tensor indices 
associated with each edge are separately summed. Such TNSs may be represented pictorially 
where we draw a "leg" for each index of each tensor and join contracted indices with lines. 
Physical indices are drawn as unpaired arrows. For example, the simple tensor network 
resulting from the multiplication of two matrices A^^S^^, is represented by: 



A 



B 



(2.3) ^ _ _ 

In the case of an MPS we associate with each tuple of matrices A^'', regarded as a three-index 
tensor the diagram according to 

(2.4) 



3k 



A 



In the pictorial representation the coefficient of \i1i2 ■ ■ ■ jn) for an MPS is depicted as 
(2.5) 



{ujl\A^^A^^ ■■■A^ur) 



n 



3n 



A 



A 



A 



The contraction involved in the definition of a tensor network state may also be expressed 
in terms of a path integral. To do this we define the following discrete "action" 

(2-6) S[{ai, aa . . . , a\E\); (ji, • • • ,j\v\)] = ^-i ^og{A^^^^^^^...^^^J, 

With this definition, the TNS is given by 

(2.7) \^)= I PaPje^^["'^']|j), 
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where / VotDj denotes here a discrete sum over all paths ex = (ai, 02, . . . , o:\e\) and (ji, . . . ,j\v\) 
with ak e {1,. . . ,D} and jk G {1, 2, . . . , d}. 

Intuitively speaking, the way to take a continuum hmit of a TNS is to choose the tensors 
so that as the spacing between the sites goes to zero the density of particles/excitations in 
the system remains constant. In one dimension, for MPS, this is achieved [181112] by choosing 

= I + eQ 

(2.8) = 



where Q and R are arbitrary D x D matrices and e is t he la ttice spacing. We'll see in the 
next section that with this choice of A^s the path integral (2.7) reduces, in the limit e — )■ 0, to 
a standard path integral. In the sequel to this paper, we'll show that a similar recipe works 
for any sufficiently regular lattice. 



3. Path integrals and continuous matrix product states 

Continuous matrix product states are a variational class of states for one-dimensional quan- 
tum fields. We focus on the bosonic case with field annihilation and creation operators 
and 'ip\{x) such that ["^^(x), ^^^(y)] = 5{x — y). A cMPS is then defined in terms of the 
quantum field A and an auxiliary D-level quantum system B by 



(3.1) 



\x) = (wil^exp 


7') 


f K{s) (g) I + iR{s) (g) ij\{s) - iR^s) (g) ^^(s) ds 






where K is a. D x D hermitian matrix and R is D x D complex matrix, \u)l,r) are D- 
dimensional states of the auxiliary system B, V'^(s) is a bosonic field operator on the physical 
system A, \^a) is the fock vacuum, and V denotes path ordering. 

3.1. A path integral for the auxiliary system. We wish to reformulate the cMPS state 



(3.1) so that expectation values for the auxiliary system are recast as path-integral ex- 
pressions. The motivation for this is two-fold: firstly, to facilitate the passage to higher- 
dimensional cMPS states; and secondly, to make manifest the symmetries of the physical 
state in terms of symmetries of an action for the auxiliary system. Our discussion is centred 
on the case of a single bosonic field in (1+1) dimensions; the generalisation to spinor and 
vector fields follows easily, and we only comment on the modifications required. 

Write a basis for the hilbert space Tis of as | j = 0,1, . . . , D — 1}. We enlarge this 
space via second quantisation, i.e., by either introducing bosonic annihilation and creation 
operators bj and foj according to the canonical commutation relations: 

(3.2) [bj,bi]=6,,k, J = 0,1,...,D-1, 

with all other commutators vanishing, or fermionic annihilation and creation operators Cj and 
Cj according to the canonical anticommutation relations 

(3.3) {cj,cl} = ^j,k, j = 0,l,...,D -1, 
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with all other anticommutators vanishing. The configuration space for our enlarged auxil- 
iary system is that of fock space 5^±(7^b), where the ± subscript indicates either bosonic or 
fermionic fock space. 

The connection between T-Lq and our enlarged fock space §^±("^6) is made, in the bosonic 
case, by identifying with the single-particle sector via 

(3.4) \j)s^b]\n)is, 
or, in the fermionic case, 

(3.5) \j)^^c]\n)s, 

where is the fock vacuum. We identify, whenever clear from the context, states |t<;) G 
with their single-particle counterparts in §^±(7^6). 

Using this embedding, a cMPS (3.1) is equivalent, in the bosonic case, to 



(3.6) \x) = {uL\Uil,0)\un)\nA) = {coL\Vexp 



-i [ F{s)ds 
Jo 



where F is a one-parameter set of field operators on AB, generated by U{1, 0) and given is by 

D 

(3.7) F{s) = K^\s)h]h ® 1 + iR'\s)h]hk ® i^\{s) - tR* {s)b]bk ® ^\{s), 

j,k=i 

This equivalence of definitions follows from the fact that F{s) is particle-number conserving 
on system B (i.e. its action on B is through terms of the form 6^6^ only), and so we remain 
in the single-particle sector throughout. The fermionic version is identical except that bj 
operators are replaced with cjs. 

The parameter s can be regarded as a time coordinate for the auxiliary system. We then 
obtain a path integral by dividing [0,/] into small intervals Sq = 0, Si, S2, ■ ■ ■ sn = I with 
uniform spacing Sk+i — Sk = e, so that U (/, 0) = U{l,l — e)U (/ — e, / — 2e) ■ • ■ f/ (e, 0), and insert 
resolutions of the identity between each term. Our choice of resolution is, in the bosonic case, 
in terms of coherent states of the auxiliary system, defined as \(j)k) = exp[0fc6^ — (plbkllfle)'- 

(3.8) ^ " ^ / n^''^'^ I <pk){^k<j)k\, 

^ k 

where N = l/e. In the fermionic case we exploit fermion coherent states of the form 
\(j)k) = exp[c|.0fc — (j)lck]\^B), where 0^ are now Grassmann-valued. Apart from the use 
of anticommuting Grassmann numbers the fermionic calculation mirrors the bosonic case in 
essentially all other respects; we thus focus on the details of the bosonic calculation and write 
out the fermionic case at the end. 

After the resolution (3.8) has been placed between each term we end up with a product 
of transition amplitudes of the form (®fc0fc(s + e)|?7(s + e, s)| ®fc' (pk'{s)) ~ {^k'Pkis + e)\l — 
ieF{s) \ <^k' 0fc'(s)). We then make use of the expression 



(3.9) (®fc0fc(s + e)| 0fc'("S)) = exp 



lf2\^k{s + e)\' + m\'-2<Pl{s + e)Ms) 



k= 
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and the assumption that only smooth variations of (pkis) contribute, which allows us to 
expand the terms in the exponential and obtain, in the continuum limit e — 0, 



(3.10) 



where the path integral is over D complex fields and S is an operator-valued action given by 
(3.11) S = [ ds ('i0"^(9,0 - - i{^^R^)^\ + i(0ti?t0)^O ^ 



where we've written {(pk} as a vector (f). However, since the field operator '?/'^(s) commutes 
with iPa{s') and 'ip^is') at all other points s' the ordering over the auxiliary time variable is 
trivial and we can simply write the path integral as 



(3.12) 



\X)= /l)20exp[z^(0,0t)]|$) 



where |$) is a physical field coherent state 



(3.13) 



1$) = exp 



Hs)^\{s)-<i>*{s)^A{s) ds 



$(s) = (l)^R(j), and the complex action S is given by 

(3.14) 5(0, 0t) = f ds {i(l)^ds(p - (p^Kcp) . 



This formulation (3.12) of the one- dimensional cMPS state as a path integral is one of our 



guiding expressions for the generalisation to higher-dimensional scenarios which we describe 
in a sequel paper. The fermionic case is identical, except that is now a vector of Grassmann 
fields. Of course, both the bosonic and fermionic cases yield exactly the same physical state 
\x), since they coincide on the 1-particle sector. Notice that we've dropped the limits from 



the integrals; the expression (3.12) makes equal sense for quantum systems on [0, /] as for the 



infinite case (— oo, oo) 



3.2. Interpretation of the cMPS path integral. Let's pause to interpret (3.1) for a 
moment. What these equations are saying is that a cMPS is a superposition of coherent 
states 1$) with some weighting e*'^. The standard intuition concerning coherent states is 
that they are the "most classical" states of a quantum system due to their saturation of 



the Heisenberg uncertainty relation. Thus, (3.12) tells us that a cMPS is a superposition of 
"classical" field states centred around classical field configurations $ : M — )■ C in phase space. 
These field configurations $ themselves are scalar functions of a vector of auxiliary classical 
fields : M — ?► C. By interpreting that spatial variable s as a temporal variable one can 
understand the action 5* for these auxiliary fields as that of a (0 + l)-dimensional quantum 
field, i.e., ordinary quantum mechanics. 

One has the picture of an auxiliary system undergoing a classical trajectory of its dis- 
crete variables, however to gain information (by measurement) about a dynamically evolving 
quantum system we inevitable disturb it because of the back-action of the quantum measure- 
ment. The closest representation of the dynamics in this quantum setting is to continuously 



CONTINUUM TENSOR NETWORKS 



7 




\x) 



I 



X 



$(x) = <p\x)R(j){x) 



Figure 1. An illustration of the coherent-state path integral representation 
for a cMPS state \x)'- here a sample classical trajectory for the (in this case, 
three) auxiliary fields is depicted above. These classical trajectories are then 
combined via (f)\x)R(f){x) into a single complex scalar trajectory The field 

coherent state is then represented diagrammatically via the ket with the purple 
smeared trajectory (the idea is that a field coherent state is like a smeared-out 
classical configuration centred on the classical configuration ^{x) oc Im($(a:))). 
The formula for the resulting cMPS is then a superposition of such coherent 
states weighted by the "action" S. 



monitor the evolving auxiliary system with a sequence of infinitesimally weak measurements 
[3]. By exploiting von Neumann's prescription for quantum measurement this process is then 
understood as entangling the auxiliary system and an infinite collection of meter systems. 
The combined auxiliary system+meter collection undergoes completely positive dynamics. 
In the continuum limit the meter systems constitute a quantum field with one extra spatial 
dimension, the reduced state of the meters alone is a quantum state. The cMPS coherent 
field state is then an imprint of the discrete trajectory, and is as classical a record as possible. 
The stength and manner of this imprint is entirely contained in the particular coupling R{t). 
Each trajectory for the auxiliary system contributes a coherent field state, and the cMPS is 
simply a superposition of "classical" trajectories with the according weighting by the action 



In this section we show that an arbitrary one-dimensional quantum field state is expressible 
as a cMPS, as long as the bond dimension D is allowed to become arbitrarily large. The 
argument we present here is for the case of bosonic fock space i5^_(L^([0, Z])) on a finite 
interval [0,Z], however, it is expected to be valid in the more general case of the interval 
(— oo, oo). 



S. 



4. Completeness of the cMPS class 
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The argument is rather simple and rehes on three facts. The first is that an arbitrary 
quantum field coherent state 



(4.1) |$)=exp 



Jo 



is exactly representible as a cMPS |x(^)) with bond dimension D = 1. This follows upon 
taking K and R to be the one-dimensional matrices K{s) = and R{s) = $(s). The 
boundary vectors \ul) and \ujr) are simply taken to be equal to 1. The next fact we require 
is that the set of all field coherent states is dense in Fock space, meaning that an arbitrary 
field state G ^-(-^^([0, ^])) may be approximated arbitrarily well by an increasing linear 
combination of field coherent states: 

N 

(4.2) 5^c,|$,)^|M^). 

1=0 

The final fact we need is that a linear combination \x) = ci|xi) + C2IX2) of two cMPS and 
1X2) with bond dimensions Di and D2, respectively, is again a cMPS with bond dimension 
D = Di + D2 and parameters K = Ki Q) K2, R = Ri Q) R2, = (ci(wl,i| © C2(co'l,2|), and 

Putting these facts together allows us to deduce that 

N 

(4.3) |XAr) = $^c,|<l>,) 

1=0 

is a sequence of cMPS with bond dimensions Dn = N that tend, in the limit, to an arbitrary 
state l^') in fock space. Thus we have confirmed the completeness or expressiveness property 
of the cMPS variational class in one dimension. It is worth noting that the argument we 
present here is by no means the most economical: there are, exploiting gauge invariance, al- 
most certainly more efficient sequences of representations tending to the state |\E') using lower 
bond dimensions. Indeed, as we argue in the next section, a more economical representation 
of a physical field state is strongly suggested by the path integral representation. 

It is worth noting that in the previous subsection we showed that an arbitrary cMPS 
is a superposition of field coherent states. Here we've shown the converse: an arbitrary 
superposition of field coherent states is also a cMPS. 

5. CMPS REPRESENTATION OF PHYSICAL QUANTUM FIELD STATES VIA PATH INTEGRALS 

In this section we show how naturally occurring physical states are efficiently represented 
via cMPS. Intuitively, the argument presented here is a continuum version of that appearing 
in [2], i.e., we exchange the role of space and time and contract the (continuum) tensor 
network for the state of a dynamically evolving field along the spatial axis first, regarding it 
is a temporal axis. 

We frame our results in the bosonic setting, although the extension to the fermion case is 
straightforward. The generic situation we consider is therefore that of bosons moving in M 
with field annihilation iIj{x) and creation (x) operators obeying the canonical commutation 
relations [ip{x),ip''{y)] = S{x — y), with all other commutators vanishing. 
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The hamiltonian, in second quantisation, is taken to have the standard kinetic+potential 
energy form: 

(5.1) H = f + V + W, 
where 

(5.2) f=[ ^^iM^^rfa;, v = [ Vix)i^\x)^ix)dx, 

J ax ax J 

and with an interaction potential 

(5.3) W= I W{x -y)il)\y)'4)\x)il){x)il){y)dxdy. 



For simphcity we will concentrate here on ultralocal interactions, i.e. W{x — y)= wS{x — y), 
with w being a constant. 

We consider a field system, initialized in the physical state \f{0)) G "H^ and allow it 
to evolve under the full hamiltonian H for a time T, until it reaches the state \ip{T)) = 
e~^^^\ip{0)). Our task is to show that the final state \^(T)) admits a natural description in 
terms of the cMPS path integral representation, which can be interpreted instead as a (virtual) 
process in which some additional auxiliary system undergoes dissipative dynamics that 
couple it to the physical field system and on completion generates \ip(T)). To avoid confusion, 
in the case of the physical evolution of the field system 'H_4 we use x for the spatial coordinate, 
and t for the physical time coordinate, while for the virtual process in which the auxiliary 
system T-Lb couples to the physical field we use s for the virtual time coordinate of T-Lb 
and label subsystems of with the parameter /3. The construction that follows will roughly 
amount to reinterpreting the field variables {x, t) as (s, /3) within a quite physically distinct 
setting. 

Our first move is to reformulate the physical field evolution in terms of a path integral 
expression over coherent states. The construction proceeds, as usual, by discretising the time 
interval [0,T] into n pieces of length e = T/n and writing 

(5.4) |^(T)) = (e-^)l^(O)). 

We suppose, for simplicity, that the initial state \ f{0)) is a field coherent state. 

As in the construction of the auxiliary action, we insert a resolution of the identity, in 
terms of 1-d field coherent states |$(t)) := exp[Jdx ^{x,t)'tp^ (x) — ^*{x,t)'ip{x)]\fl), between 
each application of e~"^. Expanding up to first order, and using the key overlap equation 



(3.9) we find that an infinitesimal advance for the physical system is described by 



(5.5) ($(t)|e"''-^|$(t - e)) ^ ^.-^ ^*{=^,t)dtH=^,t)-dt^*{x,t)'S>(x,t)-ien{^*{x,t)M^,t))dx ^ 
with a hamiltonian density ^{x^t) given by 

(5.6) H{^*{x,t),^{x,t)) = \d^^{x,t)\^ + V{x)\^{x,t)\^ + w\^{x,t)\^ 
Summing over each time interval yields the expression 

(5.7) \^{T)) = j P$P$*e'^(*'**)|$(T)), 



CHRISTOPH BROCKT, JUTHO HAEGEMAN, DAVID JENNINGS, TOBIAS J. OSBORNE, AND FRANK VERSTRAETE 



being a superposition of physical coherent states described by $(a;,T) at time t = T, and 
with the action 

(5.8) ^($,$*)=/ / i^*{x,t)dt^{x,t) -ni^*ix,t),^{x,t))dxdt. 

Jo J-oo 

The lower limit of this path integral is 0) = (p{x, 0) while the upper limit is unconstrained. 

The path integral form of \(p{T)) is suggestive of how an auxiliary system should couple to 
the physical system in order to generate \ip{T)) under (virtual) dissipative dynamics. Since we 
wish the auxilary system to sweep over the length of the physical field the time parameter 
for the process s, should correspond to the physical spatial variable x. 

To capture this idea we could subdivide the auxiliary system T-Lb into harmonic oscillator 
subsystems as "Hg = ^jBl-ip, labelled by some variable /3, however since t is a continuous 
variable we shall effectively be taking the limit in which "Hg is an auxiliary complex field 
where (3 is its spatial coordinate and the auxiliary system has spatial extent from /3 = 
to /3 = T. The key point is that spatial couplings (along /3) within the hamiltonian of the 
auxiliary system can be used to simulate the physical dynamics that generates |v?(T)), as 
the auxiliary system sweeps out over the physical field, and couples to it through a natural 
interaction term. 

For the auxiliary variables we use 2o(s, /3) and 2i(s, /3), which we can combine into a single 
complex field as z = + i^i- The free hamiltonian of the auxiliary system is taken to be 

(5.9) 

J 

-i(2o(s,/3) - i2i(s,/3))9/3(2o(s,/3) +iSi(s,/3))rf/3, 
and where po and pi are the momenta conjugate to 2o and zi. 



The form of (5.7) suggests that the interaction term coupling the auxiliary and physical 
systems be taken to be the continuous measurement interaction in which the physical system 
'K_A is interpreted as continuously measuring the 'observable' 'z = zq \ i'z\. This is obtained 
as the continuum limit of the coupling 

(5.10) ifi,t(s) = 5^ b{s - je) [?(s, /3 = T) ® ^J_^ - z\s, [3 = T) ® t^,^] , 

in other words, the physical system only couples to the extreme edge of the auxiliary system 
at the (auxiliary) spatial point (5 = T. Here = and a^^ is the operator which 
annihilates a boson with wavefunction ■^X[{j-i)e,je){.^) for the physical system. 

It is now a case of checking that the composite system l-Lj^ ® T-Ljg, evolving under the full 
hamiltonian i^tot = K + Hint for auxiliary time from s = —oo to s = +oo will indeed generate 



the desired field state \(p{T)) as expressed in the path integral form (5.7). The calculation 
proceeds in a similar manner to the earlier cMPS path integral calculation evolving under 
the composite hamiltonian i^tot, however for our resolution of the identity at auxiliary time 
s we use the complete set of states {|-2(s))} given by 

(5.11) \z{s)) = \zo{s, 0), zo{s, e), ■ ■ ■ zo{s, T); zi{s, e), zi{s, 2e), ■ ■ ■ zi{s, T)), 
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which we express in the discretized setting with oscillators located at /3 = 0, e, 2e, . . . , T. A 
straightfoward calculation gives that 

(5.12) (a;z.|Pe-*/^-W'^>R)|^]) = j V^z{s, ^)V'p{s, ^) exp[iS' {po,Pi, zo, z,)]\n,,) 
where we have the action 

/oo pT 
ds / dl3 {pozo + pizi - K{p, z) - iz{s, (3 = T)ijjl^ + iz*{s, j3 = T)-iPa)- 
■oo Jo 

Consequently, by identifying z{s,(3 = T) with $(a;,T) we see that the evolved physical state 



\ip{T)) can be represented by a cMPS with free hamiltonian K given by (5.9) and interaction 



given by (5.10). 



The cMPS representation that we have constructed involves an infinite dimensional auxil- 
iary system where integration over /3 corresponds to a continuum summation over the auxil- 
iary indices; this is not unexpected since the auxiliary system faithfully simulates the entire 
dynamical history of the physical field. However, the local character of the interaction term 
implies that we can obtain \(p{T)) equally well from the coupling of a single auxiliary oscillator 
to the physical field, with the composite system now undergoing a more general completely- 
positive map (instead of a unitary interaction). Specifically, the above calculation has shown 
that |(^(T)) = {ujL\U\uR)\n), or more generally |y?(T))(y?(T)| = Tt^^^[U{u ^ \n){n\)W] for 
some operator U on the joint system and auxiliary state u, but which can now be written as 
Tr^=r[Tr^^7-[f/(ci; (g) )[/"'']] = Tti3=t[^{^i3 ® |^)(^|)] for some completely-positive map 

S defined on the physical field and oscillator at /3 = T. By truncation of the oscillator hilbert 
space, and simulation of the evolution S we may thus obtain an efficient cMPS description of 
\f(T)) in terms of a purely discrete auxiliary quantum system. 



6. Conclusions 

In this paper we have investigated the cMPS variational class of quantum field states from a 
path integral perspective. Building on the the observation that discrete MPS representations 
can be viewed as path integral sums that couple discrete auxiliary and physical 'trajectories', 
we have constructed a natural path integral representation of cMPS as a sum over coherent 
field states. This representation is physically appealing, allows us to show the completeness 
of the cMPS class, and guides the construction of higher dimensional continuum limit of 
tensor network states. The construction also allows a natural representation of a dynamically 
evolved physical field state. 
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